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Abstract. The Wigner delay time is addressed semiclassically using the Miller's S- 
matrix expressed in terms of open orbits. This leads to a very appealing expression, 
in terms of classical paths, for the energy averaged Wigner time delay in chaotic 
scattering. The same approach also puts in evidence the semiclassical incapability to 
correctly assess the time delay higher moments. This limitation suggests that the use 
of the semiclassical approximation to quantify fluctuations in scattering phenomena, 
like in mesoscopic physics, has to be considered with great caution. 
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What is the time a quantum particle spends to traverse a scattering region? 
Reported superluminal wave propagation pQ and its relation to the fascinating and 
controversial problem of the "tunneling time" El H] brought renewed interest to the 
quantum (or wave) dwell time question in general. The status of this matter is very 
nicely discussed by a recent comprehensive and pedagogical review jH] • 

The present study is devoted to examine the delay time of a particle scattered by 
a chaotic potential [HI Ej in the absence of tunneling barriers. Here the semiclassical 
approximation is employed to take a fresh look over the Wigner time delay [B] . Starting 
directly from the Miller semiclassical S matrix, a quite appealing expression for the 
average dwell time is obtained and some quite unexpected limitations to the semiclassical 
approximation are clearly revealed. 

The Wigner-Smith time delay matrix (HI E] is defined as 

N qc 

Q ab (E) = -ih^s:^ , (i) 

where the scattering matrix S, that encodes all accessible information about the 
scattering process, is taken at the energy E. The sum in runs over all N open 
asymptotic scattering channels. The Wigner time delay Tw(E) is then given as 

T W (E) = ±trQ(E). (2) 

In scattering processes large enhancements of t w (E), or long dwell times, 
correspond to narrow isolated resonances. In these situations, the Wigner time delay 
gives access to properties of individual quasi-bound states. In the opposite and also 
relevant case of overlapping resonances, information about individual states is lost. The 
fluctuations of tw(E) result from the coherent interference of many resonances. Here 
the scattering process is best characterized by suitable statistical measures of Tw(E) 
averaged over energy windows containing many resonances. In chaotic systems tw 
displays universal fluctuations that can be described by the theory of random matrices. 
A very detailed discussion on various aspects on the statistical approach for the Wigner 
time delay can be found, for instance, in [TP]. 

The semiclassical literature on the Wigner time delay is invariantly based on its 
relation to the level density [TT | fT^ | . as envisaged long time ago by Friedel [T51II1]. Hence, 
the resulting semiclassical tw depends solely on the properties of the periodic orbits 
trapped in the scattering region. The semiclassical calculated Wigner time fluctuations 
coincide with the random matrix results ^HlEIj for N ^> 1, as they should. 

It is also desirable to cast the semiclassical Tw in terms of open classical orbits that 
spend a finite time in the scattering region, closer to the spirit of a scattering problem. 
This is the task pursued here. Equation (JTJ) is directly evaluated using the Miller's 
semiclassical S'-matrix [T7j . namely 

S ab {E)= £ (3) 

/i(a«— b) 

where the classical trajectories that start at channel b and end at channel a are labelled 
by //(a <— b). Accordingly, is the reduced action (with a Maslov phase included) 
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and p^ is the classical transition probability for going from b to a following the path /i 
[0]. Throughout this paper the wide tilde indicates quantities obtained by means of the 
semiclassical approximation. In the derivation of (J3J) the absence of tunneling barriers 
between the scattering and the asymptotic regions is implicit. Furthermore, the number 
of open channels N must fulfill N ^> 1. 

For the semiclassical Wigner time delay, energy variations in classical transition 
probabilities p^ are negligible as compared to those in the actions cr^, since the latter 
are measured in units of h. Hence, 

% = ^Ew^ ei( " )/7i > ( 4 ) 

where t^ = da^/dE is the classical time the particle spends to go from b to a through 
the path //. Here the sums are unrestricted and run over all classical trajectories that 
enter and leave the scattering region. It is worth stressing that, due to the semiclassical 
approximation, fw is not manifestly real for any given value of energy E. In analogy 
to the unitarity deficit of the semiclassical S'-matrix ^H], it can be shown that the 
imaginary part of fw is of sub leading order in powers of 1/N. This spurious imaginary 
part can be easily eliminated by using, for instance, 



d 



e\ 



QT(E) = -ihJ2fe S ac (E + -)S: C (E - -) (5) 



e=0 



2' acy 2' 

instead of Q a b defined by Restricting the analysis to the leading order term in 
powers of 1/N, it is possible to insist with Q ab . In addition to the simplicity, this 
strategy has the merit of exposing some of the semiclassical limitations. 

The statistical analysis of fw and its higher moments unveils system specific, as 
well as universal features [Hj. Let us start discussing the energy averaged fw, namely 

<^> = ^E(^v^^" . (6) 

Here (• • •) indicates an energy average taken over an energy window AE where the 
classical dynamics presents little changes, nonetheless comprising many resonances. 
To compute (fw) it is justified to neglect the energy dependence of the probabilities 
Pfj, and use the diagonal approximation. The latter says that, in general, different 
orbits of a chaotic system are uncorrelated, and holds for trajectories with dwell times 
shorter than the Heisenberg time th = h/A. Here A is the mean resonance spacing. 
Fortunately, without barriers, trajectories with t exceeding th are statistically negligible 
in the semiclassical regime of N ^> 1 [12]. In the absence of system specific symmetries 
the diagonal approximation reads (exp[i(cr /t — a v )/K\) = 5^, yielding 

= jtT.^Vp^ 6 ^ = JjUpJh = t ■ ( 7 ) 

This remarkable equation expresses the multichannel energy averaged Wigner time delay 
purely in terms of classical quantities: r is the classical scattering region escape time 
obtained by averaging over all trajectories weighted by their transition probabilities. 
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The energy average eliminates the quantum interference terms in the leading order. 
Equation (JJJ) holds for perfect transmission, that is guaranteed in the semiclassical limit 
by the absence of tunneling barriers. The diagonal approximation is customarily justified 
only for chaotic systems. Hence, chaos is a key element to obtain (J7|). The validity of 
the diagonal approximation for other kinds of dynamics is unclear. Thus, despite of the 
simplicity and appeal of (J7J), a rigorous derivation of an expression of the same kind for 
integrable systems is still lacking. 

This (semi) classical result raises an important question: Is (fw) = h/(NA) 
consistent with the well-known exact identity {tw{E)) = h/(NA), that holds 
irrespective of the underlying classical dynamics [Zj? In other words, is it possible 
to prove r = h/(NA) using just geometric arguments? 

The equivalence between both relations can be shown for the special case of chaotic 
billiards connected to the scattering region by wave guides. The proof is straightforward. 
The average time between bounces T b for an ergodic billiard is 

ib = for d = 2 and r b = — — for d = 3 (8) 

Pv Av 



where v = ^2E/m is the particle velocity. In two dimensions (d = 2), A/P is the 
billiard surface to perimeter ratio. Likewise, V/A is the volume to surface ratio for 
three-dimensional (d = 3) cavities. The above relations are exact. As discussed in jT9~] . 
their rigorous derivation is long known by mathematicians. Simple physical heuristic 
arguments can also used to obtained (JBj), as shown in [2T?j . 

For ergodic billiards the relation between the scape time r and the average bounce 
time Tb is 

P A 
t = — r b for d = 2 and r = — Tb for d = 3 (9) 
W S 

where W and S are respectively the wave guide width in d = 2 and its cross section 

area in d = 3. The Weyl formula can be used twice to relate the geometry of the billiard 

to the mean resonance spacing A and the number of channels in the wave guides N, 

yielding 

r = ~~ (10) 

for both two and three dimensional systems, as previously announced. 

Let us now use the "open orbits" semiclassical approximation to calculate the 
Wigner time delay autocorrelation function, namely 

C(e) = (t w (E + e)r w {E)) E - (r w (E)) 2 . (11) 

More explicitly, we compute 

= 4 E E y/PvPvPsPwt^' (e*^-'" + V-<v>) e i(t,-u)e _ (f w (E)) 2 , (12) 

a.b ^,iy(a' — b) 
"A /i',K'(c<— d) 

and compare it with the result obtained using either the stochastic approach fSlEl, or 
the "closed orbits" semiclassical theory [TT] . The analysis of the transition probabilities 
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Pn gives the answer almost immediately. For chaotic systems the transition probabilities 
follow the analogue of the Hannay-Ozorio de Almeida sum rule for open systems [2JJ 122] 

E P, = Wr e ~ tlT5t (13) 

t<t^<t+8t 

where Y.Kt^t+st is the sum of all classical transition probabilities following the 
trajectories /x belonging to the a small time interval [t,t + St], where St is classically 
small. In (j!3|) the sum has no restriction on channels. Hence, it is easy to verify that 
by replacing the sum by a time integral, the classical normalization (flux conservation) 
condition [U] 

N 

E E Pm = i , (14) 

0=1 /i(a<— 6) 

is fulfilled. Note that (jUj) must hold irrespective of the system dynamics. 

The semiclassical Wigner time delay autocorrelation function C (e) can be calculated 
using the diagonal approximation 

( e s^- CT "+V-<v)^ = s^S^, + S^S^ (15) 

and the sum rule (|13)1. First, the diagonal approximation contracts pairwise the 
summations over the orbit indices in (fT2|) . As a result (fw(E)) 2 is cancelled and only a 
double sum over the paths remains to be evaluated. By grouping the paths with similar 
traversal times we then use the sum rule (JT3~j) to transform the summations over orbits 
into time integrals. The later immediately give 

C^) = ~ 1 —-- 2 . (16) 



N 2 



i + (f) ; 



Since no special attention is payed to time-reversal symmetric paths, represents the 
semiclassical correlation function for broken time-reversal symmetry. It is noteworthy 
that although f\y{E) is not purely real, the energy average eliminates the imaginary 
part of C(e). However, (jlfij) is at odds with the random matrix result fSl ITT)] , that 
reads 

2r 2 !-(f) 2 



C(e) 



N 2 



2 > 



(17) 



in the limit of N ^> 1 and perfect transmission. The agreement is not reestablished by 
just starting with a Wigner-Smith matrix Q which in the semiclassical limit is manifestly 
Hermitian, as for instance Q sym in (JSJ). Neither the variance, nor the e dependence of 
the Wigner time correlation function defined as f^ m = (l/N)trQ sym agrees with ()17j) . 

Related problems appear when the semiclassical approximation using open orbits 
is employed to obtain the average transmission through a cavity and its autocorrelation 
functions, quantities of central interest in mesoscopic quantum coherent electronic 
transport. In this case the standard semiclassical approach fails because the S matrix is 
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not unitary. As previously discussed [TH] this problem is best characterized by noticing 
that 



Although admittedly small, the lack of unitarity is of the same order as typical 
transmission correlation functions. These spurious "unitarity" fluctuations could 
perhaps be reduced by including higher order corrections in the individual semiclassical 
S matrix elements. However, a correct description of the correlations among different 
matrix elements, neglected in the standard semiclassical approach, results to be a more 
practical way to proceed. In fact, for the specific case of transmission, the discrepancies 
between random matrix theory and the semiclassical approximation could be fixed by 
introducing proper semiclassical sum rules that impose unitarity [THj . 

The Wigner time delay problem does not seem to have such a simple solution. 
Here the shortcomings of the semiclassical approach are more severe. For instance, 
a direct inspection shows that both the variance and the functional dependence of 
{tw(E + e)fw(E)) and (fw(E + e)f^(E)) are very different. Unfortunately, even (JHJ) 
which seemed very promising to eliminate the semiclassical spurious imaginary part 
feature is of little help: Although f^ m (E) is manifestly real, r^ m (E)f^y m \E + e) is not. 

In distinction to the approach presented in this study, the "closed orbits" 
semiclassical approximation to the Wigner time delay is very successful to describe its 
fluctuations |TT|. I12j. The reason is simple. The trace over Q allows one to express the 
Wigner time delay as a density of states, see for instance [15J. This is a more convenient 
starting point for the semiclassical approximation: the density of states is manifestly 
real and since Q was already traced the unitarity problems of S do not appear. As a 
result one keeps only the contributions of closed orbits trapped in the scattering region, 
loosing all information about the channels. The small price to pay is that the simple 
physical interpretation for the average fw is missed. On the other hand, since numerous 
observables of interest in scattering problems do not involve traces over all channels, a 
more general solution is desirable. 

It is important to mention that for some classes of cross-section correlations, the 
relation between random-matrix theory j2H| and semiclassical results is in very close 
relation to the above discussion. As nicely addressed by Ref . j23j , apparent discrepancies 
between the two approaches are successfully remedied by avoiding taking singular 
semiclassical limits. Unfortunately we do not see how to adapt such ideas to our 
problem: The scattering problem considered in [21] allows for the use of a "closed 
orbits" -like semiclassical approximation that is free of the unitarity problems discussed 
here. 

In summary, our study shows that the semiclassical S matrix leads to a very 
appealing expression, in terms of classical paths, for the energy averaged Wigner time 
delay in chaotic scattering. On the other hand, it also puts in clear evidence the 
semiclassical incapability to correctly assess the time delay higher moments. This 
limitation can be attributed to the spurious imaginary part of fw and to the lack of 




0. 
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unitarity of the S matrix. The unitarity problem is quite severe and not only specific to 
the object studied in this work. Hence, before this problem is circumvented, applications 
of the semiclassical approximation to quantify fluctuations in scattering phenomena, like 
in mesoscopic physics, have to be considered with great caution. 

Acknowledgments 

The authors thank CNPq and PRONEX (Brazil) for partial financial support. 
References 

[1] Wang L J, Kuzmich A, Dogaiu A 2000 Nature 406, 277 

Steinberg A M Kwiat P G Chiao R Y 1993 Phys. Rev. Lett. 71 708 
Biittiker M and Landauer R 1982 Phys. Rev. Lett. 49 1739 
Biittiker M 1983 Phys. Rev. B 2 6178 
Nusscnzveig H M 2000 Phys. Rev. A 62 042107 
de Carvalho C A A and Nussenzveig H M 2002 Phys. Rep. 364 83 

Smilansky U 1991 Les Houches 1989 Session LII on Chaos and Quantum Physics ed M J Gianonni, 

A Voros and J Zinn- Justin (Amsterdam: North-Holland) pp 371-441 
Lewenkopf C H and Wcidcnmiiller H A 1991 Ann. Phys. (N. Y.) 212 53 
Wigner E P 1955 Phys. Rev. 98 145 
Smith F T 1960 Phys. Rev. 118 349 

Fyodorov Y V and Sommcrs H J 1997 J. Math. Phys. 38 1918 
Eckhardt B 1993 Chaos 3 613 

Vallejos R O, Ozorio dc Almeida A M and Lewenkopf C H 1998 J. Phys. A: Math. Gen. 31 4885 
Friedel J 1952 Philos. Mag. 43 153 
Balian R and Bloch C 1974 Ann. Phys. (N. Y.) 85 514 
Lehmann N, Savin D, Sokolov V and Sommcrs H-J 1995 Physica 86D 572 
Lehmann N, Savin D, Sokolov V and Sommers H-J 1995 Nucl. Phys. A582 223 
Miller W H 1975 Advances in Chemical Physics ed K. P. Lawley (New York: Wiley) Vol. 30 p 77 
Vallejos R O and Lewenkopf C H 2001 J. Phys. A: Math. Gen. 34 2713 
Chernov N 1997 J. Stat. Phys. 88 1 
Jarzynski C 1993 Phys. Rev. E 48 4340 

Kadanoff L P and Tang C 1984 Proc. Natl. Acad. Sci. USA 81 1276 
Hannay J H and Ozorio de Almeida A M 1984 J. Phys. A: Math. Gen. 17 3429 
Fyodorov Y V and Alhassid Y 1998 Phys. Rev. A 58 R3375 
Alhassid Y and Fyodorov Y V 1998 J. Phys. Chem. 102 9577 
[24] Eckhardt B, Fishman S and Varga I 2000 Phys. Rev. E 62 7867 
Eckhardt B, Varga I and Pollner P 2001 Physica E 9 535 



[2. 
[3; 
[4; 

[5; 

[e; 

[7; 

[»; 
[9; 
[10; 

[ir 
[12; 
[13; 

[14 
[15; 
[16" 
[17; 

[is; 
[19; 

[20" 
[2i; 
[22; 
23" 



